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Abstract 

We develop the method analyzing particle number non-conserving phenomena with non- 
equilibrium quantum fieldtheory . In this study, we consider a CP violating model with interaction 
Hamiltonian that breaks particle number conservation. To derive the quantum Boltzmann equa- 
tion for the particle number, we solve Schwinger-Dyson equation, which are obtained from two 
particle irreducible closed-time-path (2PI CTP) effective action. In this calculation , we show the 
contribution from interaction Hamiltonian to the time evolution of expectation value of particle 
number. 
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I. INTRODUCTION 



There are several particle number violating phenomena in the elementary particle physics 
such as leptogenesis. We are intereted in them when they occur in the non-equilibrium 
state. To understand how particle number is produced, we study a simple model with 
particle number violating interaction. We treat the interaction perturbatively and show 
how the interaction Hamiltonian contribute to the expectation value of the particle number. 
Throught out our paper, we use the 2PI CTP formalism. The 2PI CTP formalism is a 
method which leads to the Green functions for equilibrium and non-equilibrium states. The 
stationaly condition of 2PI effective action determines the Schwinger- Dyson equation for 
Green functions. The Green functions depend on the initial density matrix which represents 
the initial statistical state of the system. In the next section, we introduce our model and 
show the process which leads to the particle number violation. 

II. MODEL 

We consider the model that creates the particle number. In figure 1, we show how 
particle number can be produced in our model. The top pannel in the figure corresponds to 
the particle number decreasing process while the bottom figure corresponds to the particle 
number creating process. The difference of the rate asymmetry calculated from both of the 
diagrams corresponds to the net particle number. In the figure 1, A is a coupling constant 
which violates the particle number, while A^ corresponds to a coupling constant which 
conserves the particle number . B 2 is a mass term which violates the particle number. The 
Lagrangian which describes four processes in figure 1 is given by, 

£ = ^Nd^N-lm%N 2 + \d^\ 2 -ml\<p\ 2 +B 2 <p 2 +B* 2 <p* 2 +AN<p 2 +A (1) 

where N is a neutral scalar field and the particle number is assigned to be zero. is a 
complex scalar field which particle number is assigned to be 1. The Lagrangian contains 
two sources of particle number violation. One is an interaction term A^, and the other is 
the mass term B 2 . A CP violating phase of the present model corresponds to the phase 
difference between coupling constants A and B. It is given as, 

A = arg[A]. (2 ) 
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We rewrite the complex scalar field by introducing the two real scalar fields, <$>\ and 2 - 
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FIG. 1: Feynman diagrams which are related to the particle number violation. 
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Then the Lagrangian of eq. (OQ) is written as, 

1 1 1 2 1 2 

i=l i=l 

where mi and m 2 are masses of the real scalar fields 0i and 02- They are given as , 

mi = (mj - S 2 ) 1/2 ,m 2 = (m! + B 2 ) 1 ' 2 . 



(3) 
(4) 



(5) 
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Note that the masses of the real scalar fields are splitted into mi and m 2 by the influence 
of B 2 term. denotes the interaction of N f> 4>i(f)j and is given as, 



|A|cos[<M + ^ -\A\MM 
-|A|sin[<^] |A|cos[0 A ]-^ 

In the following sections, we use the Lagrangian in eq.fjSJ). 



(7) 



III. RELATION BETWEEN U(l) CURRENT AND GREEN FUNCTIONS 

The particle number which we describe in the previous section corresponds to the U(l) 
charge associated with the phase transformation of the complex scalar field (p. In this section, 
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we relate the U(l) current to a Keldysh Green function In general, there are four types 
of Keldysh Green functions. Their definitions are, 

G}}(x.,y)=<T < f H {x)<l> j {y)> ( 8 ) 

G%{x,y)=<f<f> i {x)<t> j (y)> (9) 

Gl](x,y)=<$(y)<j ) l(x)> (10) 

G%(x,y)=<<Pl(x)$(y)>, (11) 

where the matrix elements implies the trace with the initial density matrix p as, < A >= 
Tr[Ap\. Using Keldysh Green function, the U(l) current can be written as, 

<j,(X) >= [d xli G%(x,y) -d yfl G{l(x,y)]\ x=y=x . (12) 

In this way, we can compute the U(l) current with the Keldysh Green functions. 



IV. KELDYSH GREEN FUNCTIONS FROM 2PI EFFECTIVE ACTION 

In this section, we derive the Schwinger- Dyson equation for the interacting field theory in 
non-equilibrium enviroment. There is the 2PI CTP formalism which introdued by E.Calzetta 
and B.L.Hu [2] to treat the non-equilibrium state. This method applies to our model. The 
2PI effective action of our model is given as 

r = ^logdet^)- 1 + ~ £ [logdet^,)- 1 ] + j ^Jf^ Gfffoy) + \ £ —£j^G$(x,y) 

id — 1 I'd — 1 

+r 2 + const, (13) 

where the action S can be written, 

S[<p a ] = j d 4 x^c ab d^N a d^N b -^c ab N a N b + ^c ab d^d^^ (14) 

where C\\ = 1,C22 = —1 and the other components are 0. By taking variation of the 2PI 
effective action T, we obtain the Schwinger-Dyson equation of Keldysh Green functions. 



4 




FIG. 2: 2PI diagram included in of the 2PI effective action of eq,(|13p 
They are written as, 

i^>5 m J{x-y) = -p^(^+<)c^(a ; ,y)-y , d 1 «ES(^^(^y) 

+ / *zK%(x,z)G%(z,v), (15) 
iS 5a 8 mn 5(x - y) = -c a Jd 2 y +m 2 n )G 5 ±(x,y) - [ d 4 zGZ(x, z)J$(z,y) 



+ J d'zGZ(x,z)Ki a (z,y), (16) 

where K is a function which is associated with the initial distribution and £ is a self-energy. 
At initial time, we assume that the expectaion value of U(l) current vanishes. Then we can 
focus on the time variation of the current due to the effect of interaction Hamiltonian. The 
divergence of the expected value of the current is, 

d, < f{Xn >= - / [/i 2 2 (*°, X\ k) - eH(X°, X°, k)] , (17) 

where the f and e are the interaction terms in Schwiger-Dyson equation .Substituting the 
explict forms of them, the divergence of the current is given as, 

d,<f(Xn > = -SAvAvMJ T^/t^/* dz°G%(z ,X ,p)G%(z ,X°,-p-k) 
x (Gll(z ,X°,k)-Gll(z ,X°,k))). (18) 

When B is 0, then Gu = Gii- Therefore the current divergence vanishes. When we compute 
the right hand side of eq. (1181) in one-loop approximation, we can use the Keldysh Green 
functions which are derived from eq. (fT3l) without T<i- After the time(2; ) integration, the 
current divergence is written, 

%<nxn>=-M 1 M i J J ^ (pVj(k ^p + iow , e.^^w^"^ 

(19) 



where a; is a energy of each particle and iVy denotes a combonation of bose distributin 
functions. 

N ++ = (njv + l)(rij + l)(rii + l)- n N rijni, (20) 

N — = (rij + l){rn + l)n N - rijrii{n N + 1), (21) 

N-+ = (n N + l^-fa + 1) - n N (n 3 + l)n u (22) 

N + - = (tin + l)(^i + - n N rij(ni + 1), (23) 

where n is bose distribution function of each particle. Each combibation corresponds to the 
process shown in figure 3. I is a function which is given as, 

ti n \ i /in , / \ _ \ COs[{co>Ar(k + p) ± UjAk) ± CUj(p)}X°] - 1 

/ WJV k + p ± w<(k) ± wy(p),X°) = ^ v ^ v ,s,f\ J • 24 

Wjv(k + p) ± Wi(k) ± C^-(p) 

From eq. (fT9l) . one can easily find non-zero current divergence can be obtained under the 
presence of the CP violating phase 4>a- In eq. (119ft . the current divergence is given by the 
sum of the four diagrams shown in the figure 3 and their four inverse processes. 

< N / <t>i 
<fr N 

y> N < N 



FIG. 3: Feymann diagrams which contribute to the current divergence. Their inverse processes are 
not shown. 



V. SUMMARY 

In this paper, we compute the current divergence of our model by using 2PI CTP 
formalism. As a result, we show the mechanism how particle number is produced. We 
find that the production rate of particle number is propotional to the CP violating phase 
<pA by computing the current divergence eq. ( 1T9l) with the Keldysh Green function. The 
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current divergence is given as in eq. (I19|) and it is the sum of the contributions in figure 3. For 
numerical caluculation of the production rate, we must carry out the momentum integration 
of eq.flHD 
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